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We study the numerical simulation of supersymmetric models having a local 
Nicolai map. The mapping can be regarded as a stochastic equation and its 
\ numerical integration provides an algorithm for the simulation of the origi- 

nal model. In this paper, the method is discussed in details and applied to 
examples in + 1 and 1 + 1 dimensions. 
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I. INTRODUCTION 



A deep property of supersymmetric field theories is the existence of the Nicolai map |JTJ] 
that is a non linear transformation of the bosonic field <fi such that: (i) the transformed 
action describes a gaussian field £ with unit covariance, (ii) the Jacobian determinant of the 
transformation exactly cancels the fermion determinant. 

In some cases, the map is local: the new field £ is expressed by a polynomial in <p and 
its derivatives. When this happens, it is convenient to regard £ as a random field and the 
transformation £ = £ ((f)) clS cL stochastic equation |2|. 

A + 1 dimensional example is Supersymmetric Quantum Mechanics || which we shall 
adopt as a toy model in later discussions. In 1 + 1 dimensions, a class of models with 
local Nicolai map is that of N = 2 Wess-Zumino (WZ) models. Other examples are 1 + 3 
dimensional N = 1 super Yang-Mills in the light cone gauge and d dimensional lattice 
linear gauge theory ||. 

The existence of a local Nicolai map plays a major role in the lattice formulation of 
such models. Indeed, its discretization provides a recipe for the construction of a lattice 
field theory which retains most of the continuum symmetries. In particular, doubling of 
the fermions is automatically implemented due to relations between the boson and fermion 
propagators [0,0- Such lattice models can be studied by standard simulation techniques, 
but it is interesting to see if the underlying discrete stochastic structure can be useful for 
the purpose of numerical computations in a more direct way. 

If one succeeds in solving the stochastic equation £ = £(0) then uncorrelated config- 
urations may be obtained by generating £ samples. Moreover, explicit fermion fields can 
be avoided because their correlation functions may be expressed by means of the so-called 
stochastic identities in terms of statistical correlations of the solution 0(£) and the 
random field £. 

The actual implementation of this programme faces three difficulties: (i) specific bound- 
ary conditions must be imposed to preserve supersymmetry making highly non trivial to 
determine 0(£) (ii) there may be more than one solution (iii) the stochastic equation may 
admit no thermal equilibrium [IU]9(]. In other words, when regarded as an evolution equation 



for 0, unbounded solutions may appear with possible instability problems in the numerical 
simulation. 

In this paper, we discuss the above three problems in + 1 and 1 + 1 dimensional specific 
examples. We study general properties of the continuum and discretized Nicolai map and 
perform explicit numerical simulations to check the feasibility of the method. 

II. STOCHASTIC EQUATIONS AND SUPERSYMMETRY 

Let <f> a (t), t G [0,P], be a time dependent field obeying the equation 

Ut) = ^Mt) + w a (<p), w a (<t>) = ^-, (2.i) 

at d(p a 
where £ a is a gaussian white noise 

(UW)) = tap6{t-t'), (2.2) 
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and W(4>) is an arbitrary potential function. For the moment we do not specify the boundary 
conditions in (|2.1|) which make the problem well posed. The field index a may include spatial 



variables which we assume to vary in a finite volume. Equation ( |2.1| ) may be regarded as a 
functional change of variable <fi — > £ which can be inverted in a certain range of £ giving rise 
to many branches (j) = </>^(£), n — 1, . . . , N{^). 

As shown in H, the field model with periodic boundary conditions and classical action 



S 



is iV = 2 supersymmetric and supertraces can be computed as stochastic averages 



Str 



Q(<f))e 



-/3H 



Tr 



Q.{<j>)e 



-(3H, 



(2.3) 



(2.4) 



V<j)Vil)Vil) e- s n((f)) = J2f V£e- s{0 n((j) {n) (0) sgndet 



where F is the fermion number, H is the hamiltonian and S(£) = | j/^ £, 2 dt 



Equation ( |2.4|) may be exploited if one is able to solve numerically Eq. fl2.1|) with periodic 
boundary conditions 0(0) = 0(/3). 



III. THE NICOLAI MAP IN THE CONTINUUM 



A. + 1 Dimensions, Supersymmetric Quantum Mechanics 

In the case of SUSY Quantum Mechanics M , Equation ( ^.1[ ) is simply 

g = /(?) + £, 

9(0) = q(P), 
and the Hamiltonian H appearing in (|2.4| ) is 



H=\p 2 + \f(q) 2 ~ \m** 



dq 



0"3 



1 
-1 



(3.1) 



(3.2) 



Some properties of ( |3.1D do not depend on £ being a random field. For this reason, we begin 
by regarding (|3.1|) as an inner map in the space C^°([0, f3]) of periodic smooth functions in 
[O,0\. 

The Jacobian of Nicolai map can be computed explicitly ]11|] and it is 



c-sinh (\jytf(q(t)) 



(3.3) 



The number N(£) of solutions <fi for a given £ is an integer valued function of £ with possible 
jumps across the critical manifold 



M c = Lit) e C-([0,/3]) | j* f'(q(t))dt = j . 



(3.4) 



The geometry of M. c provides informations on iV(£); for instance, suppose that 
Cp°([0, f3])/M. c is simply connected, then iV(£) must be constant. To evaluate it, we choose 
£ constant and deduce from 

? = /(?)+£ =► (?) 2 = /(?)? + e? =► [\q) 2 dt = 0, (3.5) 
that g is also constant and, therefore, given by the roots of the equation 

f(q) +£ = 0. (3.6) 

A particularly interesting case is that of /' > 0, which implies M c = 0. In this case, in the 
open problem 

J^ftfc) = /(?(*,»)) (3-7) 
9(0, go) = 9o, 

we have 

<9g(£, g ) 



exp f f'(q(z,q ))dz> 1, (3.8) 

JO 



<9?o 

so that A(t, g ) = q(t, go) — go is a monotone function of go and we conclude that if a solution 
exists then it must be unique. 

To see this arguments working in particular cases, let us consider what happens when 
/(g) = /j,q n with fj, > and n — 1, 2, 3. 

n = 1 

This is the simplest case. Indeed, q((3, go) is a linear function of go and the problem 

g = /ig + £, (3.9) 
q(0) = q(P), 

has the unique periodic solution 

q(t) = e^g(0) + e"* C e'^^dz, (3.10) 
Jo 

9(0) = f e~^(t)dt. 

n = 2 

In this case, C£°([0, /3])/J\A c is n °t simply connected and, moreover, we cannot solve 
explicitly the problem 

g = /xg 2 + £, (3.11) 
g(0) = g(/3). 

Motivated by what happens when g and £ are constant we guess that, in general, there 
cannot be solutions for all £ and that, when there are solutions, they come in pairs. This 
follows also from the constraint on £ 
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rf 3 rP „ 

/ £(t)dt = -fi / q 2 {t)dt < 0, (3.12) 
Jo Jo 

which excludes some £ and from the fact that Eq. (|3.11| ) is a Riccati equation. If q%(t) is 
one particular solution for a given then another solution is 92 (t) where 

g 2 (*) = + 4^y, (3.13a) 

ii> + 2gi(t)w = -1, w{0) = w{P). (3.13b) 



In terms of 



F(t) = exp ^2 jf* gi(a)da) , (3.14) 



the function tu satisfies 
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,^(F(t)u;(t)) = -F(t), (3.15) 



and the periodicity condition w(/3) = w(0) gives 

W(0) = ~ F(0)-1 So F(a)da - (3 ' 16) 

Hence, a second solution g 2 can always be found except when F(/3) = 1 namely on the 
critical manifold 

P 

q(t)dt = 0, (3.17) 

of the map. 

n = 3 

This is a case where f'(q) = 3fiq 2 is always positive. The difference A(/3, qo) is a mono- 
tonically increasing function and it is easy to prove that 

lim A(/3, g ) = ±00, (3.18) 

qo—>±oo 

(Strictly speaking this holds only at the discrete level where blowing solutions do not appear). 
From this remarks it follows that there is a unique periodic solution for each periodic 

In the above analysis, £ was assumed to be smooth. When £ is a white noise it gives 
fluctuations 5q around the £ = solutions, namely q = q* with f(q*) = 0. The size of 5q 
depends on f'(q*) and (3. Due to periodicity, we expect it to approach a constant when 
(3 — > 00 and diverge as (3 — > with no regards to the stability of the fixed point q = q*. 

An illustrative solvable example is that of f(q) = \iq. From 

«(°) = e -;i - 1 \[ ( 3 - 19 ) 

we obtain 
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and indeed we find 



<<7(0) 2 > 



<?(0) 2 >^° ' 



1 

2/1 (1 — e~^) 2 '- 



m 2 ) 



2\ !3—>+oo 



2\n\ 



^3.20) 



(3.21) 



Here q* = and the parameter /i is just f'(q*) namely what we can call the tree level mass. 
In a numerical simulation, in the (3 — > oo limit, we expect to have one solution for each 
fixed point with large deviations depressed exponentially by potential barriers controlled by 
parameters like /x. In these regimes (assuming f'(q*) ^ 0) we can compute the sign of the 
Jacobian det<9£/<9g at each q = q* and use it in a neighbourhood of that point. 

B. 1 + 1 Dimensions, WZ Models 



In 1 + 1 dimensions, let z = X\ + ix 2 and <f>(z, z) be a complex field. Equation ( |2.1| ) takes 
the form 



(3.22) 



where 77 = 771 +ir]2 with 771 and 772 real independent white noises and /(0) = w(0) +iv(4>) an 
arbitrary holomorphic function of 0. The reason for the peculiar structure of ( 3.22 ) is that it 
guarantees that the associated field model (the so-called WZ model) is Lorentz covariant || . 
The explicit form of Q3.22 ) is 



di4>i + d 2 4>2 = u((j)i, 02 ) + r)u 

#102 - ^201 = -v(<Pl, 02) + ?72- 



(3.23a) 
(3.23b) 



The numerical integration of (|3.23|) is difficult because even in the simplest cases the asso- 
ciated random flow does not admit an equilibrium distribution. To see this, it is enough 
to consider a homogeneous (independent on 22) solution of (|3.23| ) without noise. It must 
satisfy 



dxi ^ ^ 



(3.24) 



If F(<f>) is a primitive of /(0) then the quantity H = InxF(0) is a constant of motion since 

dH 



dx\ \ dx\ 



H/l =0. 



(3.25) 



If we consider a function f(z) with asymptotic power behaviour ~ z n we see that the level 
curves of H are not closed and equilibrium cannot be reached. 

As in the + 1 dimensional case, for a constant 77 we have (by periodicity) 



2^ = /(0) + t? => 2 



<90 



dz 



m d-z +V d-z 



dx\ \ dx 2 
Jo 



dz 



0. 



(3.26) 
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which implies that is a constant obtained from the equation 



f{4>) +V = 0. (3.27) 

In particular, when 77 — f 0, tends to one of the zeroes of f{4>). As in the + 1 dimensional 
case, the Nicolai map is not singular at = 0* with f(4>*) = and f'(<p*) 7^ 0. Actually, an 
infinitesimal periodic zero mode A such that <fi = <fi* + A is a solution of (|3.22| ) must satisfy 

f)\ 

2 d~z = f(0 * )A ' (3 ' 28) 

and therefore 

(Add - |f (0*)| 2 )A = 0, J d Xl dx 2 (|VA| 2 + |f (0*)| 2 |A| 2 ) = 0, (3.29) 

giving A = 0. 



IV. THE DISCRETE NICOLAI MAP 

When time is discretized the interval [0, j3] is divided into N sub intervals and the Nicolai 
map is a change of variables 

(Co, • • • , 6v-i) <-> (<?o, • • • , Qn-i)- (4.1) 

where q n and £ n are the values of q and £ at times t n = n/3/N. The next step is the 
choice of a definite discretization scheme of ( |2.1| ). Then, one has (i) a discrete stochastic 
equation whose numerical properties must be studied and (ii) a precise fermion action to be 
analyzed from the point of view of the fermion doubling problem. In this Section we address 
these problems as well as the use of the so-called stochastic identities to compute fermion 
correlation functions without introducing Grassmann fields at all. 



A. Choice of the Discretization Procedure 

For simplicity, let us consider a single component field q(t) in + 1 dimensions. The 
Langevin equation 

? = /(?) + £, (4-2) 

(amt'))=D(q(t)w-t% 

where we allowed a q dependent covariance, may be discretized according to 

q n +i = q n + e [af(q n ) + (1 - a)f(q n+1 )) + y/eD 1/2 (aq n + (1 - a)q n+1 )£ n , (4.3) 

where the constant parameter a takes into account the ambiguity in the evalutation of / 
and D. We recall that a = 1 and a = 1/2 are usually referred to as the Ito and Stratonovich 
discretization schemes ||12|| . For non constant D(q), the value of a is important. Actually, 



in the e — > limit we can replace ( |4.3|) by 
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g„+i = q n + e [/(?*) + (1 - a)L>'(g n )e^] + V~eD^ 2 (q n ) £ n + 
and read the associated Fokker-Planck equation 



<9P _ 1 <9 



l-a 



dq 



(D a P) 



d_ 

dq 



(4.4) 



(4.5) 



In our case, however, D = 1 is a constant and the same continuum limit is obtained for each 
value of a. However, the equivalence of the two discretizations is not obvious at first sight 
because the Jacobian of the change of variables (44.) depends on a. Actually, the Jacobian 
of the transformation {q} — > {£} expressed by ( |4.3|) with periodic boundary conditions 
6v = Co 5 Qn — Qo is proportional to 



J = det 



/ -l-aef' Q l-(l-a)ef[ 

-1 - aef{ 1 - (1 - a)ef 2 








\ 



V 1 - (1 - a)ef' Q 











-1 - ae/^_ 1 / 



which we can compute 

N-l N-l 

j = n (-i - + c-i)^" 1 n (i 

fc=0 fc=0 

On the other hand, in the discrete gaussian action 



/n = /(?n), 

(4.6) 
(4.7) 



(4.8) 



1 



= o" X! {(?n+l - g«.) 2 + e 2 /n} - ~~ 9n)( a /n + ( X - °)/n+l) 

We need some care to take the limit e — > in the last term. However, if we write 

X(<?n+1 - Qn)(af n + (1 - «)/n+l) = 

n 

X(^+l ~ Qn) Q(/n + /n+l) + g " (/n+1 - /n)) , 

then the correct continuum limit can be read: it depends on a and is 

1 



(4.9) 



f(q)dq + 



a f{q{t))dt. 



(4.10) 



Putting this result together with J we find a final expression which is independent from a 
and reproduces precisely Eq. (|3.^). 

In numerical simulations, the most convenient choice of a is a = 1 because in this case 
may be evaluated directly from q n without solving any equation. 
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B. Doubling Problem 



As is well known, naive discrete fermion actions are plagued by the doubling problem |TB 
In a formulation based on a discrete Nicolai map the fermion action is constrained by the 
bosonic one (another manifestation of supersymmetry) and should be free from doublers. 
In this section we briefly explain how the discrete Nicolai map accomplishes this task by 
implementing naturally Wilson type lattice fermions. To see this, let us adopt the standard 
notation 



fni 



V~f n = fn-f, 



n— 1 j 



V f n — fn+1 

v 5 /„ = ^(v+ + v-)/ n , V A f n = \{V + - v-)f ri . 



The doubling problem is related to the kinetic term only and in this Section we set / = 0. 
We begin with the + 1 dimensional case: the Ito discretization of equation q = £, is based 
on the replacement q — > V + g. This fixes the fermion matrix d^/dq and leads to the fermion 
action 



(4.11) 



On the other hand, since V + = V 5 + V A , we can identify in Eq. ( |4.11| ) the naive fermion 
action plus a Wilson term 



1 i 

Sf = - ^2lp n (lpn+l ~ 1pn-l) + T^^^ni^n+l ~ 2^„ + V'n-l)- 
^ n ^ n 



(4.12) 



A similar mechanisms happens in 1 + 1 dimensions. For instance, the discrete Nicolai map 
proposed in M for the WZ models is a cubic symmetric version of q —>■ V + g. The continuum 
map 



d\(j) + ia 2 d 2 (j) = £, 
may be discretized as [] 



6 



C>2 



-i 

1 



(V+ + V^)0 1 + Vf0 2 = 6, 



(vr 



'291 



which, using = V s ± V A , leads to the fermion action 
S F = Y J ^n [Vf + ia 2 Vf + a 3 {V? + V 



1>, 



(4.13) 



(4.14a) 
(4.14b) 



(4.15) 



namely (after a redefinition of ip) the naive discrete action plus a Wilson term. 



1 With this choice of V + and V the cross terms in £| + £| cancel. 
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C. Some Exact Properties of the Discrete Maps 



In Sec. [ill A| and Sec. |L11B| we discussed some analytical tools for the study of the 
existence of periodic solutions of the Langevin equation as well as for the determination of 
their number N(£). 

Some of those conclusions hold also in the discrete case. Let us begin with the + 1 
dimensional case. The Ito discretization 



q n+ i = q n + ef(q n ) + y/e£ n: 



of the open problem ( |3.7|) gives 

dq N 
dq 

Hence, in the case f'(q) > we have 



d 



na+^u))- 



(4.16) 



(4.17) 



-r- (Qn ~ qo) > 0, 
dq 



(4.18) 



and (at least for asymptotic / ~ g 2n + 1 with positive n) it is easy to prove the existence of a 
unique periodic sequence {q n } for each periodic {£ n }- 

Concerning the WZ models, we now show that also in the discrete equations there are 
not zero modes when <fi is constant. The equation for the zero mode A is (see Eq. ( [3.28] )) 



Vf + ia 2 V s 2 + a 3 (Vf + V% - m) - a x u 2 



A = 0, 



du 



0=00 



Expanding the periodic A in Fourier series 

X i (x 1 ,x 2 ) = c iMM< 

kl,k2 



i(k 1 x 1 +k2X2) 



(4.19) 



(4.20) 



we find the determinant of the operator in ( [4.19 ) 

det [Vf + ia 2 Vl + a 3 (V^ + V 2 - wi) - °\u 2 



(4.21) 



n 



9? 9? 9/ 9^1 9 ^2 

sin K\ — sin k 2 — u 2 — [Ui + 2 sin h 2 sin — 



which is not zero unless u\ = u 2 = 0. 



D. Stochastic Identities 

The stochastic identities relate fermion correlation functions to stochastic averages in- 
volving the solution q of the Langevin equation and the noise £. In this paper we shall need 
only the simplest of them. To prove it at the discrete level we write 
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2tt 

d£ d£ 



271 



J V27T V2vr 
= / dndqp \ 
\ dqpdta ] 

However, in terms of the fermion matrix 



n(?(0) = 



(4.22) 



a/3 



0§a 



this means 

( n(g)e« ) = ( Wg)(</^ 

and, in particular, 

( ?a£/3 ) = ( (J^aP ) = ( ) , 

which expresses the fermion propagator in terms of the g-£ correlation. 



(4.23) 



(4.24) 



(4.25) 



V. NUMERICAL SIMULATION 

From the very existence of a local Nicolai map and previous discussions it follows an 
algorithm for the numerical computation of supertraces. The first step is the extraction of 
the gaussian random numbers Then, let <fr n be the field obeying a discretized version 

of (|2.1|); we must find the initial condition <bn such that 



A(0 O ) 



(0o) - 00 I 



0. 



(5.1) 



In + 1 dimensions it is easy to identify all solutions of (|5.1| ) as well as the sign of the 
Jacobian determinant of the Nicolai map. For the WZ models in 1 + 1 the problem is 
harder. However, as we have seen, at least for large separation of the zeroes <fi* of f(4>) 
we can expect to have one solution for each (f>* . Therefore we can use {(f)*} as starting 
guesses and take for det d^/dq its value at (ft*. In the simplest case f(<f) = /Mf) + g<t> 2 we have 
(p* = 0, — fi/g and the above regime is obtained for large \ij g. 

If an operator 0(<p) is averaged as in Q2.4]) over the realizations of {£} we obtain an 
estimate of the supertrace Str [Oexp(— j3H)\ and when (3 — > +oo we obtain (0|(9|0) in the 
case of unbroken SUSY. In the free + 1 dimensional case f(q) = fxq (and similarly in the free 
WZ models) we can solve analitically the discrete equations and determine the correlation 
functions with the stochastic algorithm or in field theory. 

The discrete Langevin equation is 



q n+1 = uq n + Ve£ n , 



1 + efi. 



(5.2) 
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The linear system 



u -1 ••• 



\ — i o o ••• 

is easily solved. The inverse of the first factor is 



\ 


( <70 \ 




( ^ \ 






= -Ve 




/ 


\ Qn-i J 







(5.3) 



1 



N 



N-2 



U 1 \ 

UJ 2 UJ 

1 uj n ~ 1 J 



(5.4) 



and therefore the initial value is 

go = " v^^rr (^"^o + • • • + l • 6v-i) 

The two-point function is 



(5.5) 



N-l 



C k - ( q q k ) - {uN _ 1)2 g 

and after a straightforward algebra 

2e cu^ 2 



a; 



i mod N+(i+k) mod TV 



(5.6) 



Ck = cosh lr ^ y 1 log " 



iV 



(5.7) 



If we take the limit N — > oo with e = /3/iV and introduce the time variable r = ek we obtain 

1 cosh (//(r - /3/2)) 



and 



Str(ge-^ge-^) = lim ( q q k > , , 

w-^oo 2/i smh(/i/3/2) 



lim Str^e-^ge"^) = (0|g(0)g(r)|0) = — exp(-//r). 



(5.8) 



(5.9) 



From the field theoretical point of view, we are computing Str(ge rH qe f3H ) with the free 
action 



S=U dt(q 2 + f , 2 q 2 ), 
I Jo 

and periodic boundary conditions. The generating functional is 

Z = ^ exp Q drJq^j ^ = exp drdr'J(r)G(r - t')J(t')^ , 



(5.10) 



(5.11) 



where 
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1 oo 

k =-°° /J 2 + 



1 



Af_k 2 
~J 2 



exp 



2irik 



(5.12) 



and 



Str(ge- r V-^ 
If we use the summation formula 
exp(— iak) 



2G(r). 



fc=— oo 



7T 



£; 2 + A 2 "Asinh(TrA) 



cosh (A (a — 7r)) 



< a < 27r, 



(5.13) 



(5.14) 



we obtain again the result of Eq.( |5.8| ). A similar computation can be carried out for the 
fermionic propagator. From the solution of ( |5.2j ) we have 



(?fcfo>- 



k-l 



which can be combined to give 

which has the correct continuum limit 

(<z(t)£(t')> = ') 



(?o6fe) 
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N-k-l 



N 



-r oo 

- y 



1-UJ 



1 

1 - 

1 



2nin 
>° // 



\/e 1 — u; 



0(Z - fc) 



N 



(5.15) 



r - r 



/ 2nin 
exp I — — [r — r 



(5.16) 



(5.17) 











associated to the fermion propagator in + 1 dimensions and confirming the stochastic 
identity. 

In the interacting case, 0o must be determined by some iterative algorithm. In 1 + 1 
dimensions it is fundamental to start from a good guess. We use the Newton-Raphson 
algorithm WM to solve iteratively the set of non linear equations 



A; 



d N ) A°) 



0. 



(5.18) 



where <jy^ is the i-th spatial component of at the A;-th time slice. The correction 



do) 



in 



do) 



do) 



is given by 



OA 



-A,, 



(5.19) 



(5.20) 
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The scheme is made more robust by introducing a relaxation parameter ujr in the update 

# - 4° ) + lu r < co R < 1. (5.21) 

Another help against numerical instabilities is to require 

< p or simply \6<f>®\ < p, (5.22) 

where p is a minimum correction threshold. The choice of the optimal % and p must be 
done empirically, but we did not find it to be critical. 

Finally, another general trick which is useful to improve numerical stability is to follow 
a bootstrap procedure and solve the problem on a L x (T — 1) lattice to provide a guess for 
the L x T problem. 

As a numerical test of the algorithm we measure the boson and fermion propagators 
in simple + 1 and 1 + 1 dimensional cases. In 1 + 1, to gain statistics, we average over 
the spatial dimension and sum over all pairs of time slices with fixed temporal separation. 
Moreover, we explicitly symmetrize the propagators under r — > (3 — r. 

The simplest interacting system in + 1 is 

f(q) = -M-gq 3 , P,9>0. (5.23) 

where dynamical breaking of supersymmetry does not occur. In Tab. (|) we show the lightest 
mass as a function of g at p = 4, j3 = 5 on a T = 200 lattice. It is obtained by fitting the 
boson and fermion propagators where the latter is computed by means of the stochastic 
identity (|4.25|) . We also show the 0(g 2 ) perturbative value 



E 1 = p+^-g-^-n 2 + 0(g 3 ) (5.24) 



2p J 2p 3 

In Fig. [I], just to give an example at the critical point p = we plot the two propagators at 
g = 1 and p = computed with j3 = 10 on a T = 100 lattice. As expected, the slope of the 
logarithmic plots is the same. 

In 1 + 1 dimensions, we simulate the WZ model with /(</>) = p<p + g<p 2 - In Fig. 0, we 
show the boson and fermion propagators evaluated at p = 4, g = 0.1 on a 20 x 90 lattice 
with e s = 0.1 and e t = 0.01 where e s and e t are the space and time discretization steps. The 
continuous line is the fit with the Ansatz 



C(t) = A cosh 



(5.25) 



On a 20 x 50 lattice we have varied g with the results reported in Tab. (|rj) together with 
the one loop value of the doublet mass which is 

m(g)=p-^ 9 -. (5.26) 

Finite time step errors can be investigated in a first approximation by studying numerically 
the finite lattice propagator integrated over space. This, in standard notation, is 
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D(k) =J2 e Mk ttPt)^r—- 2 (5.27) 



where 

27T 2 1 

Pt = ~rn t n t = • • • T - 1 and p t = — sin -e t p t . (5.28) 

One can study at fixed /3, the difference in the propagator as e t is varied. We checked that 
at fi = 4, e t = 0.01 and T = 50 the finite step effects are negligible. 

We also remark that with this particular choice of f(<f>) we can sample a single zero {(f)* = 
0) without violating supersymmetry. The reason is that the shifted field — {"Pi ~ /VflS ^2) 
obeys the same equations as but with \x — > —fi. For (00) and the symmetrized (ipip) this 
change has no consequences. 



VI. REMARKS AND CONCLUSIONS 

In this paper we have shown that the existence of a local Nicolai map in supersymmet- 
ric models has useful consequences for numerical computations. It allows the formulation 
of a simulation algorithm which generates statistically independent field configurations by 
solving a Langevin equation with periodic boundary conditions. The so called stochastic 
identities can be exploited to avoid Grassmann fields. The method is feasible and consistent 
numerical results are obtained in + 1 dimensions and also in 1 + 1 WZ models even if with 
some constraint in parameter space. Further developments are possible in the direction of 
more robust integration schemes for the Langevin equation as well as in the application to 
more realistic models. In particular, work is in progress on cases where the Nicolai map is 
determined perturbatively [15| and on QCD^ [|],[5| where the Jacobian of the local Nicolai 
map is constant. 
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FIGURES 



FIG. 1. Boson and fermion propagators for the + 1 dimensional model with drift f(q) = — q 3 
at P = 10.0 on a T = 100 lattice. Apart from the different normalization the slopes of the two 
logarithmic plots are equal. 

FIG. 2. Boson and Fermion propagators for the WZ model with drift f(<fi) = n4> + g4> 2 & t the 
point \x = 4, g = 0.1 on a 20 x 90 lattice with space and time steps e s = 0.1, e< = 0.01. 
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TABLES 



TABLE I. Lightest boson (ms) and fermion (mj?) masses as functions of g at \i = 4, (3 = 5 on 
a T = 200 lattice in the + 1 dimensional model with drift f(q) = —fiq — gq 3 



9 rn B 


mp 


ITlpert 


0.0 4.00(3) 


4.01(3) 


4.0000 


0.1 4.04(3) 


4.04(3) 


4.0368 


0.2 4.07(3) 


4.08(3) 


4.0750 


0.3 4.11(3) 


4.11(3) 


4.1062 


0.4 4.14(3) 


4.14(3) 


4.1388 


0.5 4.17(3) 


4.18(3) 


4.1699 


0.6 4.20(3) 


4.21(3) 


4.1997 


0.7 4.23(3) 


4.24(3) 


4.2281 


0.8 4.26(3) 


4.27(3) 


4.2550 


TABLE II. Lightest boson (mg) and fermion (mp) masses as functions of g at fi 


= 4, p = 0.5 


on a 20 x 50 lattice for the WZ model associated to f{4>) 


= M0 + # 2 - 




9 rn B 


mp 




0.0 4.00(5) 


4.00(5) 


4.000 


0.1 4.00(5) 


3.99(5) 


3.999 


0.4 3.99(5) 


3.99(5) 


3.985 


0.6 3.98(5) 


3.98(5) 


3.965 


1.2 3.96(5) 


3.97(5) 


3.859 
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